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Abstract. We propose a method to derive the line-of-sight magnetic flux density from 
measurements in the chromospheric Ca II IR line at 854.2 nm. The method combines two 
well-understood techniques, the center-of-gravity and bisector methods, in a single hybrid 
technique. The technique is tested with magneto-static simulations of a flux tube. We ap- 
ply the method to observations with the Interferometric Bidimensional Spectrometer (IBIS) 
installed at the Dunn Solar Telescope of the NSO/SP to investigate the morphology of the 
lower chromosphere, with focus on the chromospheric counterparts to the underlying pho- 
tospheric magnetic flux elements. 
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1. Introduction 

Advances in high-resolution observation tech- 
niques have yielded data suggesting the 
existence of a weak- field domain below 
the classical canopy (Wo ger et al.l l2006h . 
Magnetic fields play a major role in the 
morphology of the chromosphere such as 
it is seen in the Ha line core. The re- 
gion below referred to as "flu ctosphere" 
dWedemever-Bohm & Wogerl I2008I) o r "clapo- 
tisphere" dRutten & Uitenbroek! 199 lb contains 
presumably only weak-fields. In numerical 
simulations, this domain is generated by inter- 
fering (acoustic) shock waves that are excited 
in the photos phere and propagate in to the lay- 
ers above (cf . ICarlsson & Steinll 19941) . 
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ICauzzi et all d2008l) have shown that the 
Ca II infrared triplet provides a convenient 
and accessible diagnostic of the chromosphere. 
Confirming the shock behavior predicted by 
the radiation hydrodynamic numerical models 
in regions with no or at least weak magnetic 
field, it has been found that there is a strong in- 
fluence of magnetic field on acoustic processes 
such as shock signature surpressio n within that 
regime of the solar atmo sphere dJudge et alj 
feOOlUVecchio et al]|2009l ). From their findings 
it has been concluded that this influence may 
be larger than generally expected, in particular 
because it is likely that there exist at least weak 
fields in this domain. 

For this reason, knowledge of the chromo- 
spheric magnetic field topology and dynam- 
ics is of significant importance for the under- 
standing of the solar atmosphere. Yet, these 
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properties are poorly known due to limitations 
set by current observational techniques, in par- 
ticular when gathering data of weak-field re- 
gions in the quiet Sun. In addition, the interpre- 
tation of chromospheric spectro-polarimetric 
data has proven to be difficult as many com- 
plicating effects have to be taken into account 
such as NLTE conditions. 

Here we present and interpret high spa- 
tially resolved spectro-polarimetric measure- 
ments in the fluctosphere above strong photo- 
spheric magnetic features. 

2. Model 

A disadvantage of th e robust center of gravity 
(COG) method (e.g. IUitenbroekll2003l) to de- 
rive LOS magnetic flux is the lack of height 
resolution. To recover the height distribution of 
the magnetic flux, a bisector analysis is often 
employed but it is susceptible to noise. Thus, to 
enhance the reliability of a topology analysis, 
we suggest a technique that combines bisector 
and COG method. 

2.1. The Hybrid-COG method 

The following steps are performed in our algo- 
rithm: 

1 . For both, the 1+ V and I- V spectra, find the 
two points (A a and Aj,) of intersection of 
the line profile with a horizontal line (the 
points at which the red and blue line wing 
have the same intensity). 

2. Compute the COG of the I+V and I-V 
spectra (using the two acquired wavelength 
ranges) with 

£*d^{a+vx^)-(i+v)o*)} 

J Aa b dA {(I+V)(/t fl ) - a+YKXft 

for I+V and analogously A- for I-V. 

3. Compute the value for the magnetic flux 
density along the line-of-sight: calculate 
the difference of the COG positions of the 
I+V and I-V profiles, using 
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where Aq is the central wavelength of the 
line, A± are the wavelength positions of 
the centroids of I+V, gi is the effective 
Lande factor, and e and m are the elec- 
tron char ge and mass, res pectively, in SI 
units (see Uitenbroek 2003, and references 
therein). 

Step 1 allows to derive the magnetic flux den- 
sity at different formation height ranges, as it 
is done in the genuine bisector approach. The 
even cuts are necessary to avoid bias in the 
measurement of the COG value. Obviously, 
the technique will work best for spectral lines 
without an emission peak, as it would make 
Step 1 (and Step 2) difficult by creating four in- 
stead of two points of intersection and also bias 
the computation of the COG. Asymmetries in 
the Stokes V profiles, such as those found e.g. 
in a detaile d analysis of a Quiet S un region near 
network by Piet arila et alJ (120071) on a common 
basis, pose a further problem. They can lead 
to biased measurements of magnetic flux. To 
test the proposed method, we have applied it to 
synthetic profiles that inhibit similar behavior. 

2.2. Check of method 

A magneto hydrostatic model of a funnel-like 
expanding flux tube served as test case for the 
feasibility of our approach. The foot point of 
the flux tube in the deep photosphere has a di- 
ameter of 500 km and a magnetic field strength 
of 2500 Gauss. For the Ca II IR line, synthetic 
NLTE spectra for all four Stokes parameters 
were computed. These served as input for the 
hybrid-COG method. The resulting output is 
compared to the initial flux tube model (Fig.Q]). 

The topology of the flux tube in photo- 
sphere and chromosphere as well as the abso- 
lute magnetic flux density was recovered prop- 
erly, and adjacent intensity levels show no sig- 
nificant jumps in the result thus exhibiting lit- 
tle noise. This is likely due to the inclusion of 
many points in the computation of the COG 
as opposed to the genuine bisector method. 
However, the three-dimensional topology in 
the higher layers cannot be directly interpreted 
because a bisector does not necessarily corre- 
spond to a thin atmospheric layer. 
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Fig. 1. Result of the accuracy test of the proposed hybrid bisector-COG method, a) input model, 
and b) recovered flux density. The height scale in b) was derived using calculations of = 1 . 
The gray scale encodes the magnetic flux density and is the same for both panels. 



To be able to compare our results with the 
model input, we computed the effective forma- 
tion height for the Ca II infrared line by aver- 
aging the geometrical heights of the locations 
where = 1 over each of the wavelength in- 
tervals computed in Step 1 . 

The line formation height decreases from 
line core towards line wing and continuum. 
Furthermore, the formation height is affected 
by the presence of magnetic fields: due to the 
"Wilson effect" = 1 is different in mag- 
netic and non-magnetic regions. As an exam- 
ple, compared to the ambient plasma the opac- 
ity within the model flux tube is lowered by the 
strong magnetic field. This leads to the forma- 
tion of T| = 1 at a height that is about 200 km 
lower in the flux tube than in the non-magnetic 
atmosphere, which has been compensated for 
in Fig. [[]b). 

Comparing both panels of Fig. Q] allows 
to test the reliability and estimate the error of 
the proposed technique. Line properties such 
as formation height, and input parameters to 
technique such as the wavelength interval for 
Step 1 limit the recoverable height interval. 
The height distribution of the magnetic field 
of the modeled flux tube was recovered quali- 
tatively. However, regions below the "canopy" 
of the flux tube show a weak signal in the re- 
constructed topology with strengths of up to 
(6+1) G. Because of discontinuities in the mag- 
netic field along the LOS, the Stokes V spec- 
tra show spikes and thus bias the I+V and I-V 
spectra. This leads to a fake magnetic signal in 
the reconstructed topology. 



3. Data 

In this section, we present observations of a 
quiet Sun region located at disk-center us- 
ing the Call infrared line at 854.2 nm on 
September 22, 2008. Several persistent mag- 
netic bright points are visible in the field of 
view, forming a network element (Fig. [2] a)- 
b)). 

In dual-beam spectro-polarimetric mode, 
the IBIS narrowband channel was config- 
ured to scan the Ca II IR line at 854.2 nm 
using six modulation states with exposures 
of 90 ms. An overall cadence of about 1 
minute was achieved by sampling the line with 
30 wavelength step that were separated by 
4.3 pm, where the full width at half maximum 
(FWHM) transmission of IBIS at this wave- 
length is 4.6 pm. The output data was a data 
cube with two spatial and one spectral dimen- 
sion. Unfortunately, even after averaging 3 full 
scans after data reduction, the signal level of 
the Stokes Q and U component was below the 
noise level - no linear polarization was mea- 
sured. 

To gather context data about the magnetic 
field in the photosphere, a spectral region that 
includes both the Fe lines at 630.1 and 630.2 
as well as the telluric blend separating these 
two lines was observed with a temporal offset 
of 1 minute. We have scanned the line with 30 
wavelength steps with a step width of 3.2 pm, 
where the FWHM of IBIS at this wavelength 
is 2.3 pm. 
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3.1. Data Calibration 

For the broadband images, the standard cal- 
ibration method of dark- and gain-table cor- 
rection was applied. In a subsequent step, 
the images were reconstructed with a speckle 
interferometric code adapted for use with 
high-order adaptiv e optics corrected data 
(Wo ger etalJ l2008). The reconstructed images 
were used in a later step to reduce atmospheric 
differential image motion in both IBIS' broad- 
and narrowband channels. 

As detailed in ICavallinil ((2006), IBIS' nar- 
rowband channel consists of two Fabry-Perot 
interferometers, which are located in a col- 
limated beam. This optical setup produces 
data that requires an elaborate procedure to 
take into account the blueshift in e ach pixel 
of a single narrowband exposure dCavallinil 
2006). This procedure has been detailed by 
Janssen & Cauzzil d2006). 

The calibration of the polarimetric data 
is accomplished with procedures originating 
from those of the Advanc e d Sto kes Polarimeter 
(ASP) (e.g. iLites et all Il993h adapted for 
the IBIS instrument. Residual crosstalk be- 
tween the Stokes parameters has been re- 
moved manually. The quality of IBIS spectro- 
polarimetric data has recently been compared 
with that of the spectrograph (SOT/SP) on- 
board the HINODE satellite The result is that 
both instruments deliver very similar spectro- 
polarimetr ic data, which indic ates accurate cal- 
ibration dJudge & et al.ll2009h . 

The adapted procedure to calibrate IBIS' 
spectroscopic and spectro-polarimetric data 
will be des cribed in detail in a forthcoming 
publication jTritschleJl2009h . 

3.2. Application of Hybrid-COG 

Analysis of our data set with the suggested hy- 
brid COG method has been accomplished us- 
ing 101 intensity values ranging from 18.5% 
(Ca II IR line core) to 32% of the local con- 
tinuum intensity. To reduce salt and pepper 
noise introduced by the measured spectra, a 
median filter with 5x5 pixels (corresponding 
to 0.85x0.85 arcs 2 ) is applied to the data. 



The COG method was applied to the Fe I 
line data to gather information about the un- 
derlying photospheric magnetic flux density 
which amounted to about 350 Gauss within the 
strongest feature. The small-scale magnetic el- 
ements were likely not resolved. The Fe I line 
is formed within a relatively thin atmospheric 
layer on the Sun, rendering the COG method 
sufficient, whereas the Ca II IR line has contri- 
butions from layers covering 1000 km. 

The thus computed map of magnetic flux 
density is also a check for the values retrieved 
from the COG method using the Ca II IR line, 
which delivers a value of 250 Gauss using a cut 
at ~50% of continuum intensity. As expected, 
the formation height at that level appears to be 
higher than that of the Fe I line. 

Using the suggested hybrid COG method, 
starting at a level of 32% of the local contin- 
uum intensity, this value drops to 125 Gauss. 
This result is not surprising because a restric- 
tion of the method to the line core is equiva- 
lent to a restriction to higher layers. Further re- 
sults of our hybrid bisector-COG analysis can 
be viewed in Fig.[2]d)-f). 

Two situations are of interest in the field 
of view. At location A in Fig. [2] d), a patch 
of likely unresolved LOS magnetic flux disap- 
pears as when including only the few points in 
the line core in the hybrid COG method, and 
thus restricting the analysis to the highest lay- 
ers. It appears as if the magnetic flux does not 
reach the corresponding heights. On the other 
hand, at the location marked B in Fig. [2] d), a 
diffuse patch starts to show a filament structure 
at the high layers. 

Overall, it appears that (i) the strong pho- 
tospheric magnetic structure in the FOV frag- 
ments, and (ii) the flux density along the 
LOS, which corresponds to the vertical com- 
ponent of the field (the data were observed 
at disk-center of the Sun), becomes smaller 
and weaker. A magnetic funnel expanding with 
height would show similar characteristics: the 
magnetic field becomes increasingly horizon- 
tal with height, and the vertical flux component 
decreases as the the horizontal component, that 
is unfortunately currently not measurable with 
the IBIS instrument, increases. 
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Fig. 2. a) Continuum intensity around 852 nm, b) Ca II IR (854.2 nm) line core intensity, c) LOS 
magnetogram computed from Fe I (630.2 nm). Cuts through the recovered 3D topology of the 
magnetic flux density by applying the hybrid bisector-COG method to the Ca II IR line scan. 
Panel d) is at 32% of the local continuum intensity (photospheric), e) at 20.2% and f) at 18.5% 
(chromospheric/"fluctospheric"). The subpanel in f) shows a contrast enhanced close-up of the 
fibrils marked by the box. At location A, a patch of magnetic flux along the line-of-sight seems to 
disappear with increasing height (also see subpanel). Location B is an example for the appearance 
of a filament like structure with height. 
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4. Conclusions 

The focus of this contribution is the deriva- 
tion of the absolute LOS magnetic flux from 
spectro-polarimetric observations of the Fel 
(630.2 nm) and CalllR (854.2 nm) lines. 

We suggest to compute the 3D topology 
of the absolute LOS magnetic field flux per 
pixel using a new technique that is based on 
understood methods. The technique has been 
analyzed for accuracy for observations in a 
CalllR line in a Quiet Sun region near disk 
center. We were capable to recover the topol- 
ogy of a strong magnetic structure extending 
from photosphere to low chromosphere, and 
find that the magnetic flux becomes weaker 
with height developing a filament like structure 
as its height increases. The suggested method 
is not capable to infer filling factors. In an un- 
resolved magnetic element, the field strength 
thus cannot be recovered. 

While our findings suggest the existence 
of a funnel-like structure in the chromosphere, 
the location of the horizontal magnetic field re- 
mains unclear. Overall, we believe that a static 
flux tube funnel is a model too simplified to re- 
flect realistically the conditions present in the 
chromosphere. 

Understanding the chromospheric energy 
balance requires spatially highly resolved mea- 
surements of the magnetic field in the chromo- 
sphere. Such data, in spite of advanced instru- 
mentation available at modern telescopes to- 
day, still lack the signal-to-noise ratio to gain 
the accuracy needed to detect weak fields at 
the temporal resolution which is implied by the 
variations seen in chromospheric diagnostics. 
The only way to address this problem is build- 
ing a new generation of telescopes with larger 
apertures. 

Interpretation of weak magnetic field mea- 
surements in the chromosphere is at least of 
equal difficulty as their gathering. New tools 
need to be developed to better understand the 
chromosphere as an important layer between 
the photosphere and corona; for example, the 
height origins of the Stokes signals are unclear 
and do not only depend on contribution func- 
tion but also on height distribution of the mag- 
netic field strength. Discontinuities and steep 



gradients in magnetic field strength remain to 
be fundamental problems: they are difficult to 
detect and to recover from measurements with 
methods available today. 
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